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Abstract. Context. In the framework of the interstellar polycyclic aromatic hydrocarbons (PAHs) hypothesis, far-IR skeletal 
bands are expected to be a fingerprint of single species in this class. 

Aims. We address the question of detectability of low energy PAH vibrational bands, with respect to spectral contrast and 
intensity ratio with "classical" Aromatic Infrared Bands (AIBs). 

Methods. We extend our extablished Monte-Carlo model of the photophysics of specific PAHs in astronomical environments, 
to include rotational and anharmonic band structure. The required molecular parameters were calculated in the framework of 
■ the Density Functional Theory. 

' Results. We calculate the detailed spectral profiles of three low-energy vibrational bands of neutral naphthalene, and four 

, low-energy vibrational bands of neutral anthracene. They are used to establish detectability constraints based on intensity 

ratios with "classical" AIBs. A general procedure is suggested to select promising diagnostics, and tested on available Infrared 
Space Observatory data for the Red Rectangle nebula. 
. Conclusions. The search for single, specific PAHs in the far-IR is a challenging, but promising task, especially in view of the 

forthcoming launch of the Herschel Space Observatory. 

2 ' Key words. Astrochemistry — Line: identification — Molecular processes — ISM: individual objects: Red Rectangle — ISM: 

' lines and bands — ISM: molecules 

1. Introduction In lMulas et alJ(l2006bl) we presented and validated a proce- 

' dure to model in a systematic way the photophysics of PAHs in 

^ I photon dominated regions. In that work, we computed the com- 

■ - - The presence of polycyclic aromatic hydrocar bons (PAHs) in Plete far-IR emission spectrum of a sample of 20 molecules 

the interstellar medium (ISM) was proposed by'Leger & Puget' and their cations in three radiation fields, covering some typi- 



fl984.) and Allamandola et al. 0985), to account for the cal astronomical environments in which AIBs are observed. We 

so-called "Aromatic Infrared Bands" (AIBs), a set of emission concluded that the main problem for the detection and identifi- 

bands observed near 3.3, 6.2, 7.7, 8.6, 11. 3 and 1 2.7 um. in nation of such bands is likely to be spectral confusion and poor 

many dusty environments excit ed by UV photons (Leger etal.1 contrast against a strong background continuum: perpendicu- 



Il989t lAllamandolaeni]ll989l) . Such "classical" AIBs do not l^r bands, which are expected to display sharp Q branches, are 

permit an unambiguous identification of any single PAH, since much favoured with respect to shallower parallel bands, 
they just probe specific chemical bonds and not its overall Our approach can also model molecular rotation and thus 

structure. Indeed, despite the impressive amount of work de- obtain the expect ed band rota t ional envelopes of both far-IR 

voted to this subject over the years, no t/e^«/f/ve spectral iden- emission bands (Joblin etal. 1 l2002h and visible absorption 

tification of any specific individual member in this class ex- bands ( Malloci et a l. 2003). It can also estimate the contribu- 

ists to date. On the other hand, every single such PAH ought tions of fundamental and hot bands, with their distribution of 

to show a unique spectral fingerprint in the far-IR spectral re- anharmonic shifts, which affect the expected band profile. This 

gion, which contains the low-frequency vibrational modes as- makes modelling much heavier: if J,nax is the maximum angu- 

sociated with collect ive oscillations of the whole skeletal struc- lar momentum which is significantly populated, J^ax (Jmax + 1 ) 

ture of the molecule (Zhang et al. 1996; 'Langhoff'l996 HSalamal is the number of rotational levels which must be traced, dur- 

ll999:.JobUn et al..2002:.Mulas et al..2003>.2006a.b) . ing the simulation, for each vibrational mode. Since for the 
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molecule s considered Jmax is a few hundred jRouan et al J 19971 
lMulasll9 98). this makes the simulation Jmax (Jmax +1) ~ 10 
times more demanding from a computational point of view. 
Furthermore, it requires the knowledge of the effective rota- 
tional constants as a function of vibrational state and of the 
anharmonic correction for hot bands. 

To elucidate the role of spectral shape on band detectabil- 
ity, we present here, as a test case, a detailed calculation for 
three low-frequency bands of the simplest PAH, namely neu- 
tral naphthalene (CioHg), and for four low-frequency bands of 
the next larger PAH in the group of oligoacenes, i. e. neutral an- 
thracene. (C14H10). In particular, for naphthalene we selected 
the lowest energy a-type, b-type and c-type transition, whose 
fundamental wavelengths are calculated to be respectively at 
15.82, 27.74 and 58.56 yum; for anthracene, we selected the 
lowest energy a-type (43.68 jim) and c-type (21.26, 26.35 and 
1 10.23 jum) bands. 

In the following Sect.|2]we present our modelling approach. 
Results are then presented in Sect.|3land discussed in Sect.|3 

2. Modelling approach 

Our Monte-C arlo modelliri^ procedure is de s cribed in de- 
tail elsewhere ( Mulas" 1998: Jo blin et aDl2002t iMalloci et alJ 
l2003i. .Mulas et al.. 2003. 2006a.bl) . as well as the basic molec- 
ular parameters it requires. We here use exactly the same 
photo-absorption cross-sections and vibrational analy ses of 
neutral napthtalene and neutral anthracene as in Mulas et alJ 
where their applicability and accuracy is discussed 
in detail. Since this is a proof-of-concept work, we restricted 
ourselves to one single well-known environment, i. e. the Red 
Rectangle (RR) halo, u sing the same radiati on field (RF) pre- 
viously adopted for it in Mulas et al. ('2006j|), derived from the 
observational work of Viih et al.. (.2005.) . 

We here additionally studied in detail the distribution of an- 
harmonic shifts arising from the superposition of a large num- 
ber of hot bands on top of the fundamental band, in the stochas- 
tic process of PAH relaxation via IR emission. Moreover, we 
also calculated detailed rotational profiles obtained from statis- 
tical equilibrium. 

To do this we needed the previously unavailable anhar- 
monic correction to vibrational terms and rotational con- 
stants, and the Coriolis vibration-rotation coupling, which 
we c omputed using the Gaussian03 quantum chemistry pack- 
age ("Frisch et al.' 2003). This implements the Van Vle ck per- 
turb at ive approach to the above problem (e. g. .Clabo et alJ 
Il988l) . which we applied in the framework of the Density 
Functional Theory, with the exchange-correlation functional 
B3LYP (Becke 1993: .Stephens et al. ,1994) and the 4-31G 
gaussian basis set llFrishetalJll984 . While basis set conver- 
gence is not achieved yet at this level of theory, it is good 
enough to yield sensible results. Indeed, benchmark calcula- 
tions showed that smaller basis sets are required to obtain good 
perturbative corrections than are needed for equilibrium ge- 
ometry prop erties and harmonic vibrational analyses cBarone. 
I2004l2005l) . 

The anharmonic vibrational analysis yields the sec- 
ond-order anharmonic corrections to the vibrational energy 



levels, which are represented, in wavenumbers, by the formula 
(see e. g.jBarone 2005, and references therein): 

Evib(n 1 , . . . , Ha,) = ^0 +^ {rii + ^ j '^^i^'^ ("' ^) ("^ ^) ' 

where xo and y, , are t he vibrational a nharmonic constants as 
defined i nlBaron e & Min ichinollll995l) with the corrections re- 
ported in Barone (2004). The values of Xij which we obtained 
are listed in Table s \K1] and ! A . 31 in the Appendix for the specific 
bands we considered here. 

Accidental near degeneracies, when they occur, lead to ad- 
ditional corrections (i. e. Fermi and Darling-Dennison reso- 
nances), which are not systematic and must thus be evaluated, 
if relevant, on a case by case basis. For the three specific bench- 
mark bands which we modelled here for neutral naphthalene, 
the levels from which emission is estimated to occur signifi- 
cantly (fundamental and first few hot bands) are unperturbed by 
such resonances. For bands for which they do occur, in the case 
of neutral naphthalene and anthracene they are smaller than an- 
harmonic corrections anyway, so that the qualitative conclu- 
sions we may draw from our test cases remain valid. 

To leading order, selection rules are the same as in the har- 
monic approximation, i. e. they connect levels in which a single 
IR-active mode changes by one quantum. The energy diff'er- 
ence 

AEk (ni, . . . , Hat) = Eyib (...,«<: -H 1, ...) - Evib (. . . , n^, . . .) 

for a transition in which only the ^* quantum number changes 
asn^H-l nit is hence given by 

AEi: («i, . . . , «iv) = (Zi^ H- 2xkk (nk + 1) + ^^'-^ ' 
which may be rewritten as 

AE<: (ni, . . . , Hat) = AEk (0, . . . , 0) -H 2xkknk + '^Xitni- 

This equation shows that anharmonicity produces a shift in the 
hot bands, with respect to the fundamental, which is propor- 
tional to the vibrational quantum numbers in the lower energy 
state involved in the transition. Since IR emission by PAHs re- 
sults from a vibrational cascade, a given band will result in the 
superposition of the fundamental plus a possibly large num- 
ber of hot bands, their contributions depending on the detailed 
statistics of the process. This in turn results in a distribution of 
anharmonic shifts, which convolves the pure (non-thermal) ro- 
tational envelope of the band. Both effects are sampled by our 
Monte-Carlo procedure. 

On top of this, each single transition, due to higher or- 
der perturbative effects, turns into the superposition of transi- 
tions including a distribution of nearby states. This produces 
a lorentzian envelope, i. e. lifetime broadening, with a width 
which roughly scales with the density of vibrational states at 
the energy at which the transition occurs. Since bands in the 
low-energy modes are emitted near the end of the vibrational 
cascades (see Figs. ^ and below or slightly above the de- 
coupling energy (we assumed Ejec - 0.32 eV for neutral 
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Fig. 1. Probability distribution of the excitation energies at 
which photons are emitted in the three bands considered of neu- 
tral naphthalene in the RR halo. 

naphth alene and 0.21 eV for neutral anthracene. iMulas et alJ 
we neglect lifetime broadening in our present simu- 
lation. Specifically, two of the three bands we modelled for 
naphthalene and three of the four for anthracene are emitted 
virtually only below Edec; only the highest energy band we con- 
sidered for each molecule is emitted in a non negligible frac- 
tion at excitation energies above Edec- Moreover, comparison 
with experimental spectra of gas-phase naphthalene at room 
temperature show a rather good correspondence of the band 
near -12.8 yum with a the oretical spectrum we calculated under 
these same assumptions ("Pi rali et alj2006l) . 

As to the change of the effective rotational constants as a 
function of vibrational state, to leading order they have a linear 
dependence on vibrational quantum numbers, i. e. 

Aeif (ni, . . . , nw) = Aeff (0, . . . , 0) - ^ a,n,- 

and equivalent relations for Befj and Ces . The a,, b, and c, con- 
stants (the so-called a matrix), which include three leading 
contributions from harmonic, anharmonic and Coriolis terms in 
the molecular Hamiltonian, are almost all very small for naph- 
thalene and anthracene, producing changes in the range of a 
few parts in a thousand in the effective rotational constants. We 
report them in Tables lA!2l and lAT4l in the Appendix. 

Only in two cases, for naphthalene, the accidental near res- 
onance of the two vibrational modes respectively at 15.82 and 
15.90 yum and at 19.49 and 19.54 jum produces large Coriolis 
coupling terms between them, resulting in a ~ 4% change in 
the affected rotational constant as a function of the number 
of quanta in that vibrational state. In the case of anthracene, 
the largest change in the rotational constants is less than 0.9%, 
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Fig. 2. Same as Fig.[2for the four bands considered of neutral 
anthracene. 



again due to an accidental near resonance between the modes 
at 25.45 and 26.62 yum. 



3. Results 

3.1. Calculated band profiles 

The three bands of neutral naphthalene for which we modelled 
the detailed rovibrational and anharmonic structure are one of 
each type, i. e. an a-type, a b-type and a c-type band respec- 
tively. They are two in-plane and one out-of-plane (the latter 
being the so-called "butterfly" or "floppying") bending modes, 
and each of them is the lowest frequency transition of its type. 
For neutral anthracene we modelled the four lowest frequency 
bands, which are, in order of decreasing energy, two c-type, 
one a-type and another c-type. 

Since the lowest energy butterfly mode of both molecules 
falls in the spectral range of the Long Wavelength Spectrograph 
(LWS) on board ISO, we convolved their synthetic spectra with 
the average resolving power of LWS with (R^8250) and with- 
out (R^ 175) the Fabry-Perot filter, assuming a velocity disper- 
sion < 5 km s Similarly, since all other bands fall in the 
spectral range covered by the Short Wavelength Spectrograph 
(SWS) of ISO, we convolved their synthetic spectra with the re- 
solving power of SWS with (R^30000) and without (R^1500) 
the Fabry-Perot filter, with the same velocity dispersion. For 
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Fig. 3. Detailed rotational structure of the band at 15.82 fim of 
neutral naphthalene in the RR halo. Strong Cariolis coupling 
with the IR-inactive band at 15.90 fim causes one of the ro- 
tational constants to differ by ~ 4% between the upper and 
lower vibrational states involved in the transition, producing a 
blue shaded rotational profile. The superposition of hot bands, 
which has a red-shaded envelope, acts in the opposite direc- 
tion, so that the overall profile including both effects is almost 
symmetric. The two panels show the band as would be seen 
with ISO-SWS with and without the Fabry-Perot filter The 
vertical dotted line marks the position of the origin of the fun- 
damental band. 



details on SWS and LWS, see their manuals on the official ISO 
web page'. The resulting spectra are shown in figures |3]to|9l 

There is clearly a qualitative difference between perpendic- 
ular and parallel bands, i. e. the presence of a central Q branch 
containing about 20% of the total flux in the band. Such dif- 
ference is very apparent in the bands at 58.56 yum of naphtha- 
lene and at 26.35 yum of anthracene: in these cases, the rota- 
tional constants are very nearly the same in the two vibrational 
modes involved in the transition, which makes the resulting Q 
branches extremely sharp, with a full width at half maximum 
(PWHM) of the order of ~0.15 cm Moreover, in the absence 



http : //www . iso . vilspa .esa.es 
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Fig. 4. Same as Fig.|3]for the naphhalene band at 27.74 fim. 
The rotational constants are almost unchanged in the vibra- 
tional transition, resulting in a very symmetric rotational en- 
velope, the asymmetry being mainly due to the superposition 
of slightly displaced sovratones of the band. 

of shading^, neither the P nor the R branch is inverted, and 
the Q branch sits in the gap between them, with little contam- 
ination. The superposition of the fundamental and hot bands 
produces a well defined pattern of anharmonic shifts, clearly 
visible in the zoomed inset in the lower panel of Figs. |5] 
and 121 With the low resolution mode of LWS, the fundamen- 
tal and hot bands merge in an unique, unresolved, red-shaded 
Q branch, with a FWHM corresponding to the resolving power 
of the instrument. The P and R branches have a FWHM of the 
order of ~8 cm"'. 

In the band at 2 1 .26 of anthracene, instead, the larger differ- 
ence in rotational constants between the states involved in the 
transition and the pattern of anharmonic shifts of the hot bands 
effectively masks the central Q branch, which is just barely vis- 
ible in the high resolution panel of Fig.|6l 

- A difference in the rotational constants between the upper and 
lower states involved in a transition leads to the inversion of some 
branches, with the consequent formation of band heads, and results in 
an overall asymmetry in the rotational profile. A band is then com- 
monly said to be red-shaded when it has a more extended tail on the 
red side, blue-shaded in the opposite case tHerzberg. 1991 ajati 
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Fig. 5. Same as Fig.|3for the naphthalene band at 58.56 jum. 
The rotational constants are almost unchanged in the vibra- 
tional transition, resulting in a very symmetric rotational enve- 
lope, with a Q branch standing out very clearly. The two pan- 
els show the band as would be seen with ISO-LWS with and 
without the Fabry-Perot filter In the lower panel, the higher 
resolving power allows the resolution, in the Q branch, of the 
separate contributions of different hot bands, with slightly dif- 
ferent anharmonic shifts (zoomed inset). 
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Fig. 6. Detailed rotational structure of the band at 21.26 yum 
of neutral anthracene in the RR halo. The rotational con- 
stants are almost unchanged in the vibrational transition, re- 
sulting in a very symmetric rotational envelope, with a visi- 
ble Q branch whose width is due to the superposition of hot 
bands with different anharmonic shifts. The two panels show 
the band as would be seen with ISO-SWS with and without 
the Fabry-Perot filter The vertical dotted lines marks the posi- 
tion of the origin of the fundamental band. 



The parallel bands of naphthalene are resolved in only two 
well-defined features, which in the b-type one at 27.74 jum in- 
clude a split Q branch dHerzberdl 1 99 1 bt) . The width of these 
partially resolved features is approximately ~6 cm"'. As to the 
parallel band of anthracene at 43.68 jum, relatively strong red 
shading and anharmonicity combine to produce a single, un- 
resolved, asymmetric feature with an approximate FWHM of 
~6 cm '. 

These differences in band shapes, which depend on the 
interplay of many molecular parameters, have direct conse- 
quences on their detectability. 

3.2. Comparison with observations 

We here examine in detail the case of the RR nebula, which 
is one of the reference targets for the observation of AIBs and 
thus a prime candidate for the identification of specific PAHs 



in space. ISO data for this object are available in the online 
archive. Under optically thin conditions, the observed spectrum 
of a given band produced by a given PAH, for instance neutral 
naphthalene, is 



dF t/Fcont r , ,„ dP 
dA dA J ^ dA 



(1) 



where is the underlying continuum spectrum, is the 

power isotropically emitted in the band by one molecule for the 
assumed RF, Wnaph. is its number density and the double integral 
is along the line of sight and over the solid angle observed. 
As explained in detail in Mulas et al. ( 2006b) . the IR spectrum 
emitted by a species in a given regime scales linearly with the 
intensity of the RF. Since extinction in the RR halo is negligible 
(.Viih et al..2005.) . the RF will change only due to dilution with 
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Fig. 7. Same as Fig.|6lfor the anthracene band at 26.35 fim. The 
rotational constants are almost unchanged in the vibrational 
transition, resulting in a very symmetric rotational envelope, 
with a sharp Q branch standing clearly out in the middle. In 
the lower panel, the higher resolving power allows the resolu- 
tion, in the Q branch, of the separate contributions of different 
hot bands, with sUghtly different anharmonic shifts (zoomed 
inset). 



the position in the nebula and such scaling relations are valid 

dP 

in this case. Hence, we can factor — as 

dA 



— = A (r, Q) , 



dA 



dA 



where A (r, Q) is an adimensional scaling factor, independent 
of A, which takes into account the variation of RF intensity with 

dPref 

position in the source and is independent of r and Q. We 



dA 



may thus rewrite Eq. Q as 

J" dr dO. Wnaph. A (r, Q.) 



dF _ dF^ont dPref 

dA ~ dA ^ dA 



"I" , , A naph. 7 



dA 



dA 



(2) 



where Tnaph., defined above, has the dimensions of a column 
density times a solid angle and is independent of A. Classical 
AIBs are commonly considered as typical tracers of PAHs. 



Fig. 8. Same as Fig. |6lfor the anthracene band at 43.68 yum. 
The relatively large change in the rotational constants in the 
vibrational transition results in a visibly red-shaded band. 



Using Eq. (|3 and integrating over one AIB, we obtain 



I 



dFn, 



^JdF_u^, 
AIB \ dA dA 



— '^naph. I 

Jai: 



dA 

J AIB dA 
If naph. ^AIB 



'7AIB F\iB , 



where t/aib is the fraction of the flux in the AIB produced by 
the given molecule and !Paib is the result of our model for the 
reference RF and the chosen band. Solving the above equation 
for Tnaph., we get 



If naph. - lAlB 



Faib 
Pmb 



(3) 



Let's now focus on another IR band of the same specific R^H. 
To be detectable, it must stand above the continuum by more 
than the noise level, i. e., for a 2-cr detection: 



dF, 



peak 



dF,, 



dA 



dA 



dFn. 



dA 
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Fig. 9. Same as Fig.|6lfor the anthracene band at 110.23 /urn. 
The change in rotational constants, combined with the super- 
position of hot bands, produces an asymmetric profile, in the 
middle of which a sharp Q branch clearly stands out. In the 
lower panel, the higher resolving power allows the resolution, 
in the Q branch, of the separate contributions of different hot 
bands, with slightly different anharmonic shifts (zoomed inset). 



S being the signal to noise ratio, from which follows, using 
Eq. ^ and then Eq. Q, 



[ dA j, 



T, 



naph. 



peak 



ref 



dA 



peak 



dF,, 



dF,. 



2 

> — 
~ S 



dA 



dA 



In the equation above, we implicitly assumed both this lat- 
ter band and the reference AIB are integrated over the same 
aperture on the sky. For ISO observations of the RR, this is a 
good assumption since RR is completely contained both in the 
SWS and in the LWS entr ance apertures llBregman e t al. 1993; 
iMen'shchikov et alJl2002l) . Solving for t/aib, we finally obtain 
that the band will be detected only if 



dFa 



2 

??AIB ^ 7" 



dA 



' peak 



dA 



(4) 



Table 1. Detection limits estimated for three test bands of 
neutral naphthalene, along with the parameters used to de- 
rive them. Two entries are listed for the perpendicular band at 
58.56 yum, corresponding to the resolving power of LWS with 
and without the Fabry-Perot filter. 



Band S 



(jum) 





R 




dA 


I dA jp 


w \ 




/ W 


cm^ yum / 




\ sr fim 



^33 ^12J 



peak 



15.82 
27.74 



25 
45 



4.6 10"'* 
1.5 10"'* 



58.56 35 1.5 10" 



7.3 10"-'* > 1 > 1 
1.5 10-->* > 1 > 1 

5.4 10"-' 0.52 0.42 

6.5 lO-^' 0.44 0.34 



This derivation is completely general, despite the fact that 
we wrote it for neutral naphthalene and applies to any spe- 
cific band of any specific PAH. We can now put the num- 
bers for the three bands we studied in detail for neutral naph- 
thalene. Upon examination of the previously calculated spec- 
tra, the strongest emission bands of neutral naphthalene in 
the RR are the in-plane C-H stretch at ~ 3.3 jim and the 
out-of-plane C-H bend at ~ 12.7 fj.m; the fraction of the to- 
tal IR emission of neutral naphthalene in such bands is calcu- 
lated to be respectively 58.9% and 26.5%. From observations 
of the RR available from the online ISO database, their inte- 
grated band intensities amount to F3 3 - 7.35 10"'^ W cm"^ 
and F12.7 = 3.51 10"'^ W cm"^. From our model we have, 
for neutral naphthalene, = 2.42 10"^* W sr"' and 

^^cont 



^12.7 = 9.08 10-^' W sr-', while the values of S and 



estimated from ISO archive data and those of 



dP, 



ref 



dA 



dA 
cal- 



' peak 



culated by our model (see Figs.|3to|5} are listed in Table 
along with the resulting detection limits at ~2cr level. Since 
//AIB was defined as the fraction of the flux observed in a given 
AIB which is produced by this specific molecule, a detection 
limit 7/ AIB > 1 means that the corresponding band (P' col- 
umn in the table) is undetectable with the achieved S/N ratio. 
The perpendicular band, with its prominent central Q branch, is 
a much more sensitive (about two orders of magnitude) probe 
for neutral naphthalene, yielding the lowest detection limit of 
T] ~ 0.4 for both reference AIBs considered. Resolving power 
helps a Uttle, since separating single hot bands yields a higher 
peak flux in the strongest ones. 

None of these bands are detected in ISO archive data, there- 
fore we here obtained a direct observational upper Umit on rj for 
neutral naphthalene in the RR. 

Turning to anthracene, its strongest emission bands in 
the RR are again the in-plane C-H stretch at ~ 3.3 fim 
and the out-of-plane C-H bend at ~ 11.3; the fraction of 
the total IR emission of neutral anthracene in such bands 
is calculated to be respectively 37.5% and 17.2%. The inte- 
grated intensities in the former of these bands observed in 
the RR was given above, the one of the ~ 11.3 //m band 
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Table 2. Detection limits estimated for three test bands of neu- 
tral anthracene, along with the parameters used to derive them. 
Two entries are listed for the perpendicular band at 26.35 fim, 
corresponding to the resolving power of S WS with and without 
the Fabry-Perot filter. 



Band S 



(jim) 



dF„ 



dA 
W 
cm^ yum 



dA 



'73.3 '711.3 



peak 



w 



50 3.1 10-'* 
90 1.9 10'* 



21.26 
26.35 

43.68 75 4.9 10"" 
110.23 25 1.5 10""* 



4.2 10-" 

5.1 10-28 

7.2 IQ-^^ 
3.9 10--** 
6.0 10-» 
7.2 10--' 



> 1 

> 1 

> 1 

> 1 

> 1 

> 1 



0.44 
> 1 
0.86 
0.52 
0.33 
0.27 



is measured to be Fn.3 = 1.39 10 W cm ^. Our model 
yields, for neutral anthracene, 'P3 3 - 4.7 IQ-^^ W sr"' and 



estimated from ISO archive data and those of 



dA 
calcu- 



!Pii.3 = 2.15 10-^^ W sr-i, while the values of S and 

/peak 

lated by our model (see Figs.|6lto|9} are listed in Table|2] along 
with the resulting detection limits at ~2cr level. Classical AIBs 
are not the only observable quantity against which ratios, and 
consequently rj values, can be computed. Neutral anthracene, 
phenanthrene and pyrene were proposed a s possible carriers of 
the Blue Luminescence (BL, IViih et al.l [2004. 2005), a fluo- 
rescence phenomenon observed in the Red Rectangle and, sub- 
sequently, in several other astronomical sources. In a previous 
work, we demonstrated that phenanthrene and pyrene could be 
ruled out based on the failure to detect their predicted IR emis- 
sion spectra ( Mulas et al. 2006a). Anthracene remained, of the 
three, the only candidate compatible with available ISO obser- 
vations of the Red Rectangle. The strongest constraint appeared 
to be the undetected longest wavelength band at ~110 yum, 
which however contained a degree of uncertainty due to the un- 
known band profile. We here reproduce the comparison of the 
calculated band with available observations including the de- 
tailed modelling of band profile. This is shown in Fig.^| The 
line of reasoning leading to Eq. (0) can be identically retraced 
substituting the integrated BL flux in the place o f the integrated 
flux in a classical AIB (see M ulas et al.ll2006al for a detailed 
derivation). This yields t/bl > 1 for anthracene, i. e. it confirms 
that this band is undetectable with the available ISO database 
observations. It is however apparent (see Fig.llO> that this band 
is expected to be just slightly below the detection limit, and an 
increase in S/N of a factor of ~5 ought to reveal it, if anthracene 
indeed produces the observed BL. 

4. Discussion and conclusions 

The upper limits we derived are on the relative contribution of 
specific molecules to well defined classical AIBs, specifically 



3, 



4x10" 



2x10-'" - 



-2x10" 



-4x10"'' - 




105 110 115 

Wavelength (fim) 



120 



Fig. 10. Comparison between the estimated IR emission spec- 
trum of anthracene (C14H10) and an ISO spectrum of the RR in 
the wavelength range 100 -120 nm. Calcul ated spectra, under 
diff'erent assumptions (see lMulas et alJl^)06a . for details) are 
drawn in dashed, dash-dotted, and continuous lines; the con- 
tinuum-subtracted ISO spectrum is shown as a dotted line. The 
central position of the fundamental of the expected anthracene 
band is marked by a tick, which shows the effect of anharmonic 
shifts. 



at ~3.3, ~11.3 and ~ 12.7 /vm. They can be easily converted 
into absolute abundance limits, but this implies an assumption 
on their spatial distribution in the observed source and on the 
detailed assumed scaling of RF intensity. The terms depending 
on such assumptions cancel when using ratios of bands, making 
them more robust. 

Despite obtaining very similar numerical values for the up- 
per limits on 773 3 and 7712.7 for naphthalene, their relevance is 
diff'erent. The ~ 3.3 fj.m band is produced efliciently by small 
neutral molecules; the intensity of this band is suppressed in 
cations and rapidly decays with molecular size, as shown by 
Fig. ^2 This implies that the PAHs significantly contributing to 
the -3.3 fim band are a small subset of the whole population, 
namely only the small ones. The number of different species 
in this subset is accordingly very much smaller, thus < 40% of 
the small PAHs (i. e. 773 3) is a much stronger constraint than 
< 40% of all PAHs. This trend is much weaker for the flux 
fraction in the band at ~ 12.7 fj.m, as shown by Fig.^] This 
is a consequence of the well-known variation in the position 
of the out-of-plane C-H bend depending on specific molecu- 
lar parameters (e. g. solo, duo, trio modes etc.), hence chemical 
diversity produces a large scatter which dominates over size ef- 
fects. This means that the band at ~ 12.7 fim is not produced 
preferentially by a subset of PAHs defined by size, but instead 
by structure. 

As to anthracene, we obtained rather weak constraints for 
773.3, somewhat stronger constraints on 7711 3. However, Fig. 1131 
shows that larger molecules can be expected to contribute a 
considerable fraction of the flux emitted in this band, which 
means that ~40% is not such a stringent limit for 7711 3. As for 
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Fig. 11. Percentage of total IR flux emitted in the in-plane C-H 
stretch near ~ 3 .3 um, as computed b y our model for the sample 
of molecules in Mulas et al. as a function of molecular 

size. Neutral species are represented by diamonds, cations by 
crosses. For each molecule, we here considered all bands in the 
range 3.2 to 3.37 fim, which corresponds to the observed width 
of the 3.3 /zm band in the RR. 
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Fig. 12. Same as Fig. ^2 for the out-of-plane C-H bend at ~ 
12.7 yum. For each molecule, we here considered all bands in 
the range 12.3 to 13.1 jum, which corresponds to the observed 
width of the 12.7 jjm band in the RR. 



the band at ~ 12.7 fim, the band at ~ 11.3 fj.m is produced 
preferentially by a subset of molecules defined primarily by 
structure, more than by size. 

While we applied our procedure to neutral naphthalene and 
anthracene, it can be generalised to any PAH for which our 
model is applicable and to any astronomical source for which 
observations of classical AIBs (or some other suitable tracer of 
PAHs) are available. As already mentioned, all PAHs have but- 
terfly IR-active modes, which usually are the lowest frequency 
modes of each molecule and give rise to perpendicular bands 
(kMulas et al...2006b.) . The lowest frequency vibrational modes 



Fig. 13. Same as Fig.^Jfor the out-of-plane C-H bend at ~ 
11.3 fim. For each molecule, we here considered all bands in 
the range 10.9 to 11.7 fim, which corresponds to the observed 
width of the 11.3 fj.m band in the RR. 

are almost always well separated in energy (except for highly 
symmetric molecules), meaning that strong Coriolis perturba- 
tions are not likely to occur for them and they ought to display 
Q branches which can be very sharp, since the effective mo- 
menta of inertia can be very close between the upper and lower 
states of the transition. In turn, according to the present calcula- 
tions these Q branches ought to contain ~2Q% of the total flux 
in the band, the distribution among fundamental and hot bands 
depending on the detailed statistics of the process and, essen- 
tially, on the vibrational modes which are left in an excited state 
upon decoupling. The spacing between the fundamental and 
various hot bands depends critically on the anharmonic vibra- 
tional constants Xij^ which vary with the specific molecule and 
band considered. If a perpendicular band is emitted preferen- 
tially when the molecule has a very low excitation energy, then 
relatively few hot bands will contribute to it, their Q branches 
remaining well resolved and sharp, and the P and R branches 
blending and being somewhat bluiTed as a result. The same will 
happen if a few Xij are much larger than the others, producing 
a well-defined sequence of sharp Q branches. 

We here modelled in detail the profiles of one perpendicu- 
lar band of naphthalene and three of anthracene. Among them, 
three display sharp Q branches (Figs.|5l0and|9jl, in which fun- 
damental and hot bands can be separated with adequate resolv- 
ing power, the other one (Fig.|6} is red shaded to the point of 
partially washing out its structure. In this very limited sample, 
therefore, most perpendicular bands possess a combination of 
molecular parameters which yields the maximum spectral con- 
trast. On a large sample of molecules, it can be reasonably ex- 
pected that this will happen for some intense bands, producing 
the most favorable conditions for their detection. 

In general, the average rotational energy of the molecule is 
expected to scale roughly with the average energy of the vi- 
brational photons emitted (see e. g. iRouan et alj|l997t) . which 
implies that larger molecules ought to have a lower rotational 
energy which, together with decreasing rotational constants, 
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in turn translates to naiTower rotational profiles and better de- 
tectability. 

The mixture of PAHs in space probably contains a wide 
variety of different molecules, in principle in many different 
hydrogenatio n and ionisation states. However, Le Page et al] 
( I2OOII I2OO3I) showed that in any given environment the pop- 
ulation of a given PAH is expected to be dominated by one or 
two ionisation and hydrogenation states, essentially depending 
on its size and physical conditions. Furthermore, chemical se- 
lection effects will favour the most stable species over the oth- 
ers, so that at least some species may be abundant enough to 
exceed their detection limit on some rj aib (or some other PAH 
tracer). Moreover, such a detection limit will be pushed one or 
two orders of magnitude lower by the much higher sensitivity 
of the instruments on board the forthcoming Herschel Space 
Observatory and by the much reduced dust continuum emis- 
sion at longer wavelengths. 

A systematic study of PAH band profiles and IR flux ra- 
tios is necessary to select the most promising diagnostics and 
species for identification. The search for single, specific PAHs 
in the far-lR is a challenging, but promising task. 
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Appendix A: Computed molecular parameters 
from the vibro-rotational analyses 

In the following we report some of the molecular parameters 
obtained in the vibro-rotational analyses performed for neutral 
naphthalene and neutral anthracene at the B3LYP/4-31G level 
of theory. Tables lA.ll and IA.3I report, respectively, the vibra- 
tional anharmonic constants Xij for the three bands of neutral 
naphthalene and the four bands opf neutral anthracene con- 
sidered. The changes of the effective rotational constants as 
a function of vibrational state are reported in Tables IA.2I for 
naphthalene and lA.4l for anthracene. For consistency, we li st in 
the tab les the frequencies we previou sly obtained (Mulas et alJ 
l2006bl) using the NWChem package dStraatsma et alJl2003l) at 
the same level of theory. Frequencies obtained by Gaussian03 
(Frisch et al. 2003) are coincident within numerical errors. We 
give more significant digits than the actually expected accuracy 
to distinguish very close vibrational modes. 



Mulas, Malloci, Joblin & Toublanc: Searching neutral naphthalene and anthracene in the Red Rectangle 



11 



Table A.l. Anharmonic constants x/y (expressed in cm ') ob- 
tained at the B3LYP/4-31G level for the three vibrational 
modes of naphthalene neutral for which we modelled the de- 
tailed rotational structure. 



Fundamental 15.820 //m 27.731pm 58.196 jum 



( uvn) 




(cm"'") 




3.249 


-0.126 


0.004 


0.035 


3.250 


-0.123 


-0.018 


0.039 


3.265 


-0.198 


0.083 


0.012 


3.267 


-0.198 


0.032 


0.018 


3.280 


-0.156 


0.217 


-0.014 


3.284 


-0.130 


0.175 


-0.026 


3.286 


-0.198 


0.131 


0.002 


3.288 


-0.178 


0.123 


-0.008 


6.179 


-1.360 


-1.051 


-0.449 


6.278 


-0.406 


-1.057 


-0.306 


6.407 


-1.347 


-0.198 


-0.239 


6.632 


-1.163 


-0.519 


-0.309 


6.836 


-0.735 


0.180 


-0.077 


6.842 


-0.659 


-0.563 


-0.240 


7.146 


-0.486 


0.061 


-0.055 


7.335 


-1.522 


-0.039 


-0.332 


7.364 


-1.765 


0.820 


-0.110 


7.875 


-0.630 


0.704 


-0.110 


7.961 


-0.305 


0.038 


-0.065 


8.273 


-0.406 


-1.356 


-0.327 


8.527 


0.011 


1.270 


0.207 


8.551 


-0.344 


-0.482 


-0.010 


8.658 


0.044 


0.165 


-0.055 


8.838 


-0.382 


0.441 


-0.229 


9.812 


0.082 


0.824 


-0.034 


9.906 


0.089 


0.099 


-0.066 


10.065 


-0.189 


-0.052 


-0.021 


10.141 


-0.140 


-0.115 


-0.126 


10.389 


-0.089 


-0.004 


-0.294 


10.601 


-0.024 


-0.068 


-0.157 


10.638 


-1.126 


0.150 


-0.083 


11.294 


-0.221 


-0.088 


-0.457 


11.937 


-0.325 


-0.089 


-0.245 


12.545 


-0.088 


-1.264 


0.292 


12.691 


-0.064 


-0.176 


-0.155 


12.957 


-0.335 


0.050 


-0.190 


13.312 


0.859 


-0.252 


-0.077 


13.937 


0.046 


0.087 


-0.093 


15.825 


0.105 


-0.104 


-0.200 


15.895 


-0.226 


-0.221 


-0.463 


19.503 


0.057 


-0.095 


0.339 


19.536 


0.050 


0.254 


0.807 


20.878 


0.237 


0.226 


-0.054 


21.162 


-0.225 


-0.324 


-0.349 


25.745 


0.259 


-0.295 


-0.579 


27.737 


-0.104 


0.236 


0.363 


53.950 


-0.113 


-0.460 


-0.230 


58.565 


-0.200 


0.363 


-0.258 



Table A.2. The a,, b,- and c, constants (in units of 10 cm ') 
as obtained at the B3LYP/4-3 IG level for neutral naphthalene. 
These quantities express the change of the rotational constants 
as a function of the vibrational mode. For comparison, the rota- 
tional constants A, B and C in the vibrational ground state are 
0.103289, 0.040779 and 0.029241 cm-\ respectively. 



Fundamental 


a; b; Ci 


Oum) 


(10-" cm-i) 



3.249 

3.250 

3.265 

3.267 

3.280 

3.284 

3.286 

3.288 

6.179 

6.278 

6.407 

6.632 

6.836 

6.842 

7.146 

7.335 

7.364 

7.875 

7.961 

8.273 

8.527 

8.551 

8.658 

8.838 

9.812 

9.906 

10.065 

10.141 

10.389 

10.601 

10.638 

11.294 

11.937 

12.545 

12.691 

12.957 

13.312 

13.937 

15.825 

15.895 

19.503 

19.536 

20.878 

21.162 

25.745 

27.737 

53.950 

58.565 



0.346 
0.347 
0.320 
0.323 
0.325 
0.323 
0.313 
0.313 
0.800 
0.471 
1.061 
0.661 
0.410 
-0.124 
-0.199 
1.918 
0.683 
0.299 
0.214 
-0.055 
-0.625 
0.250 
-0.199 
-0.488 
0.809 
0.782 
0.916 
0.954 
0.448 
-2.638 
2.807 
0.546 
0.194 
-2.053 
2.656 
0.333 
0.514 
0.319 
-40.137 
40.437 
-0.220 
-0.220 
0.053 
0.692 
0.733 
-1.792 
0.339 
2.085 



0.080 
0.077 
0.087 
0.081 
0.096 
0.093 
0.087 
0.088 
0.519 
0.387 
0.307 
0.241 
0.055 
0.302 
0.136 
0.531 
0.339 
0.161 
0.094 
0.458 
0.003 
0.098 
0.049 
0.102 
-0.042 
0.009 
0.163 
0.153 
0.192 
0.159 
0.004 
0.157 
0.017 
0.193 
0.057 
0.085 
0.057 
0.107 
-0.023 
0.062 
-0.037 
0.061 
0.025 
-0.008 
-0.042 
-0.050 
0.146 
-0.262 



0.065 
0.063 
0.068 
0.064 
0.073 
0.070 
0.069 
0.068 
0.136 
0.187 
0.478 
0.303 
-2.109 
2.421 
0.081 
0.395 
0.280 
0.104 
0.220 
0.324 
0.077 
0.028 
0.069 
0.324 
0.144 
0.112 
0.022 
0.017 
0.012 
0.002 
0.141 
0.005 
-0.001 
0.117 
-0.013 
-0.029 
0.139 
-0.020 
0.113 
-0.045 
-11.605 
11.782 
-0.097 
-0.104 
-0.113 
0.069 
-0.150 
-0.249 
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Table A.3. Same as Fig. lA.ll for the four vibrational modes of neutral anthracene for which we modelled the detailed rotational 
structure. 
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3.249 


-0.079 


-0.063 


-0.027 


0.045 


9.946 


0.062 


-0.100 


0.121 
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-0.032 
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11.005 


-0.348 
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11.315 


-2.030 


-0.563 


-0.161 
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-0.011 


11.690 


-0.356 


-0.376 


-0.030 


-0.100 


6.162 


-0.331 


-0.554 


-0.338 


-0.143 


11.982 


-1.000 


-0.666 


-0.035 


-0.096 


6.183 


-0.703 


-0.614 


1.810 


-0.184 


12.563 


-0.517 


-0.324 


0.328 


0.008 


6.329 


-0.421 


-0.520 


-1.023 


-0.177 


13.043 


-3.945 


-0.560 


0.021 


0.197 


6.486 


-0.559 


-0.572 


-0.028 


-0.119 


13.173 


-0.291 


-3.714 


0.032 


0.039 


6.514 


-0.376 


-0.471 


0.081 


-0.161 


13.434 


-0.667 


0.317 


0.023 


-0.098 


6.736 


-0.137 


-0.297 


-0.128 


-0.096 


13.512 


-0.712 


-0.806 


-0.017 


-0.027 


6.860 


-0.295 


0.048 


-0.554 


-0.126 


13.709 


-0.367 


-0.668 


0.031 


0.035 


6.864 


-0.021 


-0.365 


0.243 


-0.040 


15.329 


-0.554 


-0.136 


-0.012 


-0.209 


7.147 


-0.144 


0.118 


-0.072 


-0.031 


15.706 


0.187 


0.120 


0.016 


-0.160 


7.218 


-0.192 


-0.813 


-0.453 


-0.110 


16.336 


0.066 


-0.027 


-0.481 


0.131 


7.218 


-0.829 


-0.933 


0.704 


-0.098 


17.109 


-0.112 


-0.470 


-0.277 


-0.088 


7.434 


-0.989 


-0.699 


-0.845 


-0.140 


18.710 


-0.287 


-0.037 


-0.073 


0.265 


7.605 


-0.008 


1.095 


0.019 


-0.067 


20.037 


-0.060 


-0.499 


0.050 


0.151 


7.768 


0.290 


-0.460 


0.512 


-0.031 


21.019 


0.209 


-0.171 


0.016 


-0.192 


7.848 


2.118 


-0.500 


-0.024 


-0.087 


21.263 


0.040 


-0.065 


-0.030 


-0.002 


7.899 


-0.355 


-0.226 


0.006 


-0.135 


25.448 


-0.200 


-0.338 


0.664 


-0.424 


8.318 


-0.272 


-0.491 


0.045 


-0.002 


25.625 


0.109 


-0.064 


-0.244 


0.915 


8.473 


-0.318 


-0.152 


0.056 


-0.137 


26.349 


-0.065 


0.171 


-0.271 


-0.181 


8.553 


-0.196 


0.349 


0.069 


-0.060 


37.412 


-0.237 


-0.220 


-0.408 


0.017 


8.623 


-0.535 


2.072 


-0.178 


-0.076 


43.009 


-0.129 


-0.554 


-0.978 


-0.281 


8.668 


-0.349 


-0.198 


-0.221 


-0.037 


43.681 


-0.030 


-0.271 


0.024 


-0.587 


9.107 


-0.588 


-1.000 


-1.709 


-0.117 


82.090 


-0.261 


-0.269 


-0.828 


-0.173 


9.911 


0.105 


-0.118 


0.200 


-0.040 


110.231 


-0.002 


-0.181 


-0.587 


-0.146 
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Table A.4. Same as Fig. IA.2l for neutral anthracene. The rotational constants A, B and C in the vibrational ground state are 
0.071063, 0.014959 and 0.012361 cm-', respectively. 



Fund. a, b, c, | Fund. a, b, c, 

(jim) (10-* cm-') ifim) (IQ-* cm-') 



3 


249 


0.166 


0.014 


0.014 


9.946 


0.433 


-0.001 


0.027 


3 


249 


0.166 


0.014 


0.014 


10.102 


0.430 


0.025 


0.005 


3 


263 


0.155 


0.015 


0.015 


10.116 


0.434 


0.024 


0.005 


3 


263 


0.155 


0.015 


0.014 


10.407 


0.193 


0.029 


0.003 


3 


282 


0.161 


0.020 


0.018 


10.470 


0.072 


0.023 


0.001 


3 


283 


0.157 


0.018 


0.016 


10.875 


0.154 


0.018 


0.035 


3 


284 


0.155 


0.018 


0.017 


11.003 


0.211 


0.031 


0.004 


3 


286 


0.152 


0.016 


0.015 


11.005 


-0.323 


0.009 


0.024 


3 


292 


0.176 


0.029 


0.024 


11.315 


0.401 


0.016 


0.000 


3 


295 


0.172 


0.029 


0.024 


11.690 


0.199 


0.026 


0.002 


6 


162 


0.296 


0.106 


0.069 


11.982 


0.101 


0.024 


0.001 


6 


183 


0.375 


0.147 


0.084 


12.563 


0.261 


0.052 


0.059 


6 


329 


0.206 


0.122 


0.082 


13.043 


0.141 


0.013 


-0.005 


6 


486 


0.539 


0.081 


0.106 


13.173 


0.238 


0.008 


-0.002 


6 


514 


0.356 


0.086 


0.087 


13.434 


0.223 


-0.004 


-0.010 


6 


736 


0.276 


0.055 


0.063 


13.512 


0.313 


0.016 


0.024 


6 


860 


-0.009 


0.092 


-0.239 


13.709 


0.092 


0.014 


-0.002 


6 


864 


0.198 


0.024 


0.334 


15.329 


0.008 


0.037 


-0.033 


7 


147 


-0.073 


0.036 


0.027 


15.706 


-1.289 


-0.003 


0.015 


7 


218 


0.671 


0.117 


0.117 


16.336 


-0.177 


0.011 


0.090 


7 


218 


0.784 


0.089 


0.079 


17.109 


1.376 


-0.018 


-0.015 


7 


434 


0.409 


0.136 


0.096 


18.710 


-0.250 


0.015 


0.022 


7 


605 


-0.010 


0.095 


0.087 


20.037 


0.424 


-0.019 


-0.028 


7 


768 


0.111 


0.028 


0.000 


21.019 


0.194 


-0.007 


-0.023 


7 


848 


0.082 


0.052 


0.062 


21.263 


0.013 


-0.004 


-0.023 


7 


899 


0.283 


0.136 


0.142 


25.448 


-0.683 


-0.002 


-1.019 


8 


318 


0.165 


0.009 


0.031 


25.625 


-0.012 


0.020 


1.062 


8 


473 


-0.268 


0.004 


0.014 


26.349 


0.355 


-0.023 


-0.032 


8 


553 


-0.180 


0.035 


0.037 


37412 


0.261 


-0.006 


-0.034 


8 


623 


0.328 


0.071 


-0.448 


43.009 


0.904 


-0.060 


-0.059 


8 


668 


-0.162 


0.017 


0.565 


43.681 


-1.292 


-0.026 


0.009 


9 


107 


-0.236 


0.070 


0.055 


82.090 


0.124 


0.026 


-0.042 


9 


911 


0.391 


0.002 


0.030 


110.231 


1.647 


-0.082 


-0.091 



